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Abstract 

We investigate how the bi-large mixing required by the recent neutrino data can be accommo¬ 
dated in the supersymmetric standard model allowing bilinear R-parity violation and non-universal 
soft terms. In this scheme, the tree-level contribution and the so-called Grossman-Haber one-loop 
diagrams are two major sources of the neutrino mass matrix. The relative size of these two contri¬ 
butions falls into the right range to generate the atmospheric and solar neutrino mass hierarchy. 
On the other hand, the bi-large mixing is typically obtained by a mild tuning of input parameters 
to arrange a partial cancellation among various contributions. 


PACS numbers: 14.60.Pq, 12.60.Jv, 11.30.Fs 





Recently, impressive progress has been made in atmospheric and solar neutrino experi¬ 
ments B0. They provided us convincing evidences for three active neutrino oscillations 
requiring two large and one small mixing angles | 3 j. The resulting neutrino mixing matrix 
01 takes the form; 
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where Cy = cos 0 y, s tJ = sin d tJ , and S13 ~ 6*13 < 0.2. Here we put 6*23 = 7t/4 for the nearly 
maximal atmospheric neutrino mixing angle. The solar neutrino mixing angle 0 12 takes the 
value tan 6*12 ~ 0.65 for the so-called LMA solution which is strongly favored by the recent 
SNO data 0 . The mass-squared differences explaining the atmospheric and solar neutrino 
data are A m 2 tm ~ 2.5 x 1CT 3 eV 2 and A m 2 sol «5x 1CT 5 eV 2 , respectively. Even though 
less favored, the so-called LOW solution with tan 6*12 ~ 0.77 and A m 2 sol ~ 1CT 7 eV 2 is still 
viable. 

One of attractive schemes to generate neutrino masses and mixing is to invoke R-parity 
and lepton-number violation allowed in the supersymmetric standard model ||. The pur¬ 
pose of this paper is to address the question whether the bi-large mixing of three active 
neutrinos can arise naturally from the bilinear R-parity violation. The superpotential of the 
supersymmetric standard model may contain the following bilinear terms; 

W = einLiH-2 , (2) 

generalizing the usual fi- term, hHiH 2 . Then, there are also six soft supersymmetry breaking 
terms in the scalar potential; 

Vsoft T fTT >h 1 BiHi T h.c. , ( 3 ) 

where we used the same notations for the superfields and their scalar components. Let us 
note that Bi in the first term is dimension-one and the corresponding term for the Higgs 
bilicar is [iBHiH 2 . 

If the universal boundary condition is imposed on the soft-terms, the differences between 
the soft-terms of the Higgs boson H 1 and slepton Li such as 

A Bi = B — Bi and Am 2 = m 2 Hi — m ? L ., 
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vanish at the mediation scale of supersymmetry breaking and their non-zero values are 
generated at the weak scale through renormalization group evolution (RGE), while m 2 L iHl 
remain vanishing. In this case, there are only three free parameters e* which makes the model 
very economic. However, this model cannot accommodate the bi-large mixing consistently 
with small t/ e3 . It is easy to understand it qualitatively as one can expect that the three 
parameters e* control all the mixing angles. A small 6 1 3 and a large #23 requires ei -C 62 ~ 63 
leading to 9 12 ~ $13 f6j. j 7 |. Thus, in order to accommodate the bi-large neutrino mixing, one 
has to go beyond this minimal scheme. One way is to allow trilinear couplings while keeping 
the universality. In this case, the five couplings related to the third generation fermions may 
play a major role to generate the desired neutrino mass matrix ill- Another way is to allow 
non-universal soft-terms IBB [TT| . Introduction of general flavor-mixing soft-masses is, of 
course, tightly constrained by the flavor changing neutral current processes, such as n —> ey 
or r —> fi'y [[ 1 ^]. However, the non-universality in the flavor-diagonal soft-parameters is not 
severely constrained. Generically, one could expect A mf/m 2 Hl and A Bi/B to be of order 
one. One can also have m 2 L . Hl ~ ei‘m 2 Hl . 

In this paper, we investigate how the desired neutrino mass and mixing pattern can arise 
under such a generic non-universality condition. We will see that the right values of the 
mixing angles and the mass hierarchy can be obtained in reasonable ranges of parameter 
space without severe fine-tunning. I11 the below, we will first quantify all the tree-level and 
one-loop contributions to the neutrino mass matrix and identify the dominant contributions. 
Obtaining a rather simple form of the leading neutrino mass matrix, we will make qualitative 
discussions to understand how the desired masses and mixing arise. This will be completed 
by presenting our numerical analysis. 

Let us start our main discussion by describing the structure of neutrino mass matrix 
coming from R-parity violation. Adopting the notations of Ref. |J , the most general one- 
loop renormalized neutrino mass matrix can be written as 


nru Mf 2 M| 
ML = — =^£i£jC 2 B - z 




F N ™ J p F n 
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(€i$j + fii£j) C /3 + nh , 


( 4 ) 


where F N = M 1 M 2 /Mx / + Mfc 2/ g//i with My = c^Mi + s^ v M 2 . Here, the first term is the 
neutrino mass matrix arising at tree-level, the second terms containing Si come from the 
one-loop correction to the neutrino-neutralino mixing masses projected on to the neutrino 
direction, and the last term H^ is the one-loop correction to the /y-z/, Majorana mass 
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matrix. The non-zero values of = (L°)/(//°) — e* arise due to non-universal soft terms in 
the slepton- Higgs sector as follows; 

Am? + A Biiitp rn\. Hl 
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( 5 ) 
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where the sneutrino mass-squared is m~. = + M|c 2/?/2. As is well-known, the tree-level 

mass matrix makes massive only one neutrino in the direction of £, which is typically the 
heaviest one, z/3. In fact, the quantity controls the neutrino-neutralino mixing and thus 
could be probed by lepton flavor violating decays of the lightest neutralino in the future 
colliders 0B0- Here, let us introduce another quantity, 
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( 6 ) 


which governs the mixing between the sleptons and Higgs bosons. As we will see, the flavor 
structure of the neutrino mass matrix depends on these two R-parity violating parameters, 

and r) tl as well as non-universal slepton masses. 

A simplication of the full neutrino mass matirx comes from the obseravtion that the 
second term on the right-hand side of Eq. <w can be ignored in our case 0. This can 
be seen immediately by going to the basis where the tree-level mass matrix is diagonalized 
by the eigenvector £ and any two orthogonal unit vectors. In this basis, one finds that 
the second mass matrix has vanishing components in the 1-2 plane orthogonal to £. Thus, 
leaving the heaviest u 3 untouched, approximate see-saw diagonalization can be applied to get 
the contribution to the 1-2 plane of the order of Mzd 2 . This is like a two-loop contribution 
much smaller than the (non-vanishing) 1-2 components of the last term Tl v . Thus, there is 
no need to compute the second mass term in most cases even though we included it in our 
analysis. 

The main contribution to the last term H u of Eq. ([]]) comes from the one-loop dia¬ 
grams exchanging sneutrinos/Higgs bosons and gauginos ] 15 ], [L6| in the case of generic 
non-universality under consideration. Here we present the explicit formula of this one-loop 
mass matrix which is calculated by the use of approximate see-saw rotation ||; 
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where N a b is the 4x4 neutralino diagonalization matrix, x° denotes the neutralino mass 
eigenstates, </> represents the neutral Higgs bosons (</> = h, H and A), and the loop-function 
B 0 is given by B 0 (x, y ) = —In | — In ^ + 1 with the renormalization scale Q. The effect 
of the bilinear R-parity violating terms are encoded in the coefficients 9 ^ and Z %3 which are 
given by 
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where ry is defined in Eq. (|J), a is the usual diagonalization angle of two CP even Higgs 
bosons, and F l s = (m~ — m A )(ml — ml) (ml. — m 2 H ). Recall that the angle a is defined by 
c 2a = c 20 (m 2 A - M\)/(m 2 h - m 2 H ) and s 2a = s 20 (m\ + M 2 z )/(m 2 h - m 2 H ). 

A few remarks are in order: (i) The coefficients 9^ are the linear combinations of 9f ^s 
defined in Eq. (9) of Ref. 0. They are related by the Higgs mass diagonalization. In Eq. (||), 
the quantity appears to include the effect of neutrino-neutralino mixing by e. t . This 
dependence can be easily understood if one goes to the basis where e* vanishes [§]. (ii) The 
same diagrams have been considered in Ref. [T3j using the mass-insertion method which 
must yield the equivalent results to ours. These diagrams involve two mass-insert ions which 
can be seen here as products of two induced R-parity odd v — cj) — x° vertices, 9^9 and as 
individual sneutrino vertices, Z %3 , which is R-parity even, (iii) Among various contributions 
in 9 i( j > 9j ( j > , the term proportional to can be absorbed into the tree-level mass term giving 
a negligible effect. The term proportional to ^ t rj 3 is suppressed due to the similar reason 
discussed before, but cannot be neglected completely, (iv) The term Z tl is nothing but 
the contribution due to the sneutrino-anti-sneutrino mass splitting induced by R-parity 
violation, a la Grossman-Haber ]Hj, and Z %3 with * ^ j comes from the effective sneutrino 
mixing vertices, v % — i7* — y°. (v) The terms with Z %3 are proportional to and 

thus give smaller contributions than the terms with lyrjj from 9^9j 0 in a reasonable range 
of parameters. However, they can give a sizable effect in general. 

Now, let us consider the other one-loop contributions and show that ([!]) dominates over 
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them in the case of the general non-universality. Among various contributions, we take the 
well-known diagram with squark-quark exchange to be compared with (7|). Considering the 
trilinear couplings induced from bottom quark Yukawa couplings h b such as A ' 33 = €ih b , one 
has 
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e i e j- 


(9) 


Taking the ratio of the above two contributions, one typically gets (i/( 0 ) ps 5x 10~ 6 tp (e/ 77) 2 
with 0 = 100 GeV, n = m~ b = 250 GeV. Therefore, @ can be neglected as far as tan/3 is 
not too large and e; ~ 7 ^. In the similar way, one can find that the other diagrams are also 
sub-leading to (0). In Ref. |T(J, a slight deviation of non-universality has been assumed to 
yield e /77 ~ 10 3 and thus (|9j) was considered as the main one-loop correction. I 11 fact, this 
is a typical situation in the case of universality. The importance of the contribution (J7[) in 


the case of large deviation from universality has been notified in Ref. [|T5| and its impact on 
viable neutrino mass matrices has been considered in Refs. li 

From the previous discussions, we can write down the leading contributions to the full 
mass matrix (0) as follows: 


M u - « 
*7 


M| ££V- 

szs j c (3 
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( 10 ) 


F n '” j ~ p 32tt 2 

where /“• derivable from Eqs. (0) and (§) is the function of the masses of neutralinos, 
sneutrinos and Higgs bosons and its flavor dependence comes from the non-universal slepton 
masses. 


We are ready to discuss how the desirable neutrino masses and mixing can be realized 
by the bilinear R-parity violation with generic non-universal soft masses. For this, we will 
take the following representative set of R-parity conserving parameters; 


t p = 5, m A = 300, /j = -250, M 2 = 2M i = 200, (11) 

throughout this paper. This choice gives the light and heavy neutral Higgs boson masses, 
rrih = 84 GeV and mu = 302 GeV, respectively. Other choices will not change the main 
features of our results. Concerning the R-parity violating parameters, we allow the general 
flavor dependence for the supersymmetric e* and soft Bj parameters. To make our discussion 
simpler, we will take rn 2 L i H 1 = 0 i 11 this paper. This would be a plausible choice for the 
minimal lepton flavor violation as it may arise due to some mechanism of generating the /i 
and eifi terms. 
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Now, let us start with the simplest case: (A) the “minimal” deviation from the univer¬ 
sality , that is, sleptons have a universal soft-mass: m 2 Hl 7 ^ m 2 Ll = m | 2 = m 2 3 . This was the 
scheme employed in the analysis of Refs. JDI, IlH]. In this case, the lepton flavor dependence 


in fij disappears and thus the neutrino mass matrix ([H]) takes the following simple form: 


M" px m x XiXj + myYiYj 


( 12 ) 


where m x = |^| 2 M|c^/Fat, m y ~ |//| 2 m^o(yf 2 /647r 2 . Here, x and y are nothing but the unit 
vectors in the direction of £ and fj, respectively. As analyzed in Ref. |17|, the mass matrix 
(H) has two non-vanishing eigenvalues, m 3 ~ m x and m 2 ~ m y s 2 , whose eigenvectors are 
in the directions of x and x x (x x y), respectively. Here the angle <p is dehned by c v = x ■ y. 
From these, one finds that the desired neutrino mixing matrix is obtained for 

1 1 


x^( 0 13 ,-^,-==), and 
y oc (S12, \/ 2 (l + k)ci 2 , V 2 kc u ), 


(13) 


with an arbitrary number k. The ratio of two mass eigenvalues is given by 

m 2 9 2 m^oF N 2 \rj\ 2 2 

m 3 ~ 32vr 2 M 2 ^l^l 2 ^' 

Note that one can easily obtain its right value to accommodate the atmospheric and solar 
neutrino (LMA) mass scales; namely, m 2 /m 3 ~ Am 2 ol /Aml tm ~ 0.16 putting m y 0 = 
Fat = 200 GeV, tp = 5, |t/|/|£| = 1 and s 2 = 1. Furthermore, the relation ( |T3| ) can also be 
arranged by an appropriate choice of two independent set of parameters £j and r}j. In the 
similar way, the LOW solution can also be easily accommodated. However, it remains to be 
seen how such an arrangement for £* and ry can be made in terms of the input parameters, 
Cj, ABi/B and A m 2 /m 2 H . I 11 order to answer this question, let us choose the following set 
of values; 

ii = (0.1,1,1), rji oc (V2t 3 , 1, -1) with t 3 = 0.65 

which give rise to the desired bi-large mixing of the atmospheric and solar neutrino oscil¬ 
lations. Note that the above choice corresponds to c w = 0. The normalization of 9 will be 
chosen to reproduce a right value of Am 2 ol /Am 2 trn ~2x 10” 2 . Since we will calculate the 
ratios of neutrino mass eigenvalues and mixing angles, we put, e.g., £2 = £3 = 1- In order to 
obtain the mass scale of m 3 = 0.05 eV, one can take an overall rescaling of R-parity violating 
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variables, £, rj and e, by factor of 5 x 10 -6 . We now give three examples realizing the above 
choice of G and rji as follows. 

(Al) Am 2 /m 2 Hl = 0.7: This corresponds to the snentrino mass, = 67 GeV, and gives 
the neutrino mass matrix, 

-2.12GC, - 0. IS //,//;. 

Therefore, the choice of r/,; = (1.1,1.2, —1.2) leads to the desired results as follows; 

/\ 777.^ 

—^ = 0.03, t/ e3 = 0.08, sin 2 2d atm = 0.99, sin 2 26 sol = 0.82. (15) 

Our choice of 0 = (0.1,1,1) and the above r/,; is realized by the following input parameters; 
e, : = (4.5,1.1, -9.5) and A BjB = (0.22, 0.18, 0.23). 

(A2) A m 2 /m 2 Hl = —1: It gives rise to mp i = 228 GeV and 


-2.12CA, - 0.089 //,//,. 


Taking r/j = (1.4,1.5, —1.5), we find 


A "i 2 tm 


0.018, U e 3 = 0.07, sin 2 2d atm = 0.99, sin 2 29 so i = 0.83. 


(16) 


The corresponding input parameters are e, = (—4.2,—3.4,5.9) and AB JB = 


(-0.31,-0.15,-0.42). 

(A3) A m 2 /m 2 Hl = 0.1: It leads to = 146 GeV and 


Mh = -2.12GA, - 0.0022//,//,. 


With the choice of r/, : = (9.2,10, —10), we get 




0.02, f/ e3 = 0.02, sin 2 26> atm = 0.98, sin 2 26 soi = 0.84, 


(17) 


and the input parameters; e* = (283,287,-334) and ABJB = (0.032,0.031,0.033). 

For the cases (Al) and (A2), our general parameter scan showed that the realistic neutrino 
masses and mixing can be obtained within the range of input parameters: 1 A) |e,| G 10 
and 0.1 G \ABi/B\ < 1 leading to |£|, \r)\ ~ 1. From the above samples, one can see that 
there need certain arrangements in the flavor structure of the input parameters realizing 
the required mixing angles. This would be the case in many class of models. In our case, 
the smallness of £i is arranged not by the smallness of | 1 but by a partial cancellation 



between two terms: A m\ ~ —ABifitg leading to AB\/B « 0.22 and —0.31 for (Al) and 
(A2), respectively. This pattern arises also in more general cases as we will see shortly. Since 
| ci | is not necessarily smaller than | € 2 , 3 1, it is favored to have ABi/B ~ A B 2 ^/B. Thus, 
a vanishingly small U e3 | cannot be naturally realized in our scheme. In the case (A3), the 
universality is maintained to a certain degree, As we can see, this requires |ej| \r]i\, |£,;| 
and a strong correlation for the fine-tunned values of \ABi/B\ -Cl. In fact, this is a 
characteristic property of the universality case where the small deviation of A irtf/m 2 Hl and 
ABi/B arises due to RGE of soft parameters. We excluded such cases in our analysis. 

Let us now relax the universality condition of the slepton and Higgs boson masses, which 
leads to the following form of the neutrino mass matrix; 


•H'' = coGA, + Cij'H'lj , 


(18) 


where Co = — 2.12 again with the choice of Eq. (|TTD and the flavor dependence in c tJ appears 
due to non-universal slepton masses. We first consider an interesting case where (B) the 
flavor independence is assumed for e* to see whether only non-universality in soft-parameters 
can be the source of the bi-large mixing. As an example, we take 
(Bl) 6 i = 1, A mf/m 2 Hl = (—3.3, —3.1, —4.3) and ABi/B = (—1.0, —2.6, —3.5): 

This gives us & = (—0.047,1.25,1.27), rg = (—1.05, —1.35, —2.23) and thus 



1 0.41 

0.48 

0.15 \ 

c ij ~ 

0.48 

0.57 

0.18 


\ 0.15 

0.18 

0.046 ) 


As a result, we get 


A ml 


= 0.017, U e3 = 0.14, sin 2 2 9 atm = 0.94, 


Sill 


2B,„I = 0.73. 


(19) 


atm 


Again, one needs a relation Am 2 ~ —A Bi/atp. We find that this case (B) is not particularly 
fine-tunned compared to the previous case (A) and can be a viable option. 

Finally, we consider (C) the most general case where we take arbitrary values of the 9 
input parameters, e*, AmJ/m^ and ABi/B } whose sizes are however restricted within the 
range of (0.1 — 10). In FIGs. 1 and 2, we present the scatter plot in term of Xi = m 2 L ./m? Hl 
and pi = ABi/B with i = 1 and 2, respectively, which generate the desired neutrino masses 
and mixing. FIG. 1 shows that a solution set are centered around the values of x\ and p\ 


9 





for which the cancellation in £1 happens as discussed before. Anther solution set is allowed 
around X\ = 3.4 or 0.4 for which the sneutrino mass is close to the heavy or light Higgs mass, 
respectively. In this region, the mixing elements (j^) and thus the coefficients c 8J in Eq. flT 8 |) 
become large to enhance the one-loop contribution. As a consequence, U e 3 can be arranged 
to be small without making £1 small. In FIG. 2, one sees that the points (x 2 , p 2 ) close to 
(xi,p\) are favored although those points FIG. 1 allowing the cancellation in £1 are excluded 
as can be expected. The plot in terms of (£ 3 , p%) is also very similar to FIG. 2. In FIGs. 1 
and 2 , we plotted only the points where the tree mass is three times larger than the loop 
mass. Here, let us remark that the one-loop mass can be even larger than the tree mass. 
That is, it is possible that the one-loop contribution proportional to is the main source 
of the atmospheric neutrino mass and mixing angle while the tree mass generates the solar 
neutrino mass and mixing angle. Even though such cases of the loop dominance cannot be 
neglected, there is a much larger parameter space allowed in the case of the tree dominance 
as one can expect. This can be seen in FIGs. 3 and 4 which plotted all the allowed points 
in terms of the induced variables £j which determine the tree mass matrix as in Eq. (|4|). 
These two figures show that there appears the pattern, |£i| -C (£ 2 ! ~ |^ 3 1, which gives rise to 
#13 <C 1 s 2 3 ~ C 23 ~ l/\/2 as shown in Eq. (p~3|) for the tree-dominance case. 

To conclude, we showed how naturally the realistic neutrino mass matrix can arise from bi¬ 
linear R-parity violation assuming non-universal soft-terms. When generic non-universality 
is allowed and tan/3 is not too large, the neutrino mass matrix is dominated by two con¬ 
tributions; the tree-level mass and the one-loop mass from the so-called Grossman-Haber 
diagrams arising due to the sneutrino-Higgs mixing. This was checked by our numerical 
calculation taking the full one-loop renormalized neutrino mass matrix. In this scheme, the 
loop-to-tree mass ratio falls naturally into the right range to generate the desired values 
for Amj 0 i /Am^ ra . Considering nine input parameters, e*, A Bi and Am-, we analyzed the 
parameter space accommodating two large (#12 and # 23 ) and one small (# 13 ) mixing angles. 
Typically, the smallness of # 13 is realized by a cancellation between the terms contributing 
to £ 1 . This was shown by some examples and also by the scatter plot of FIG. 1. Such an 
arrangement would not be a severe fine-tunning of input parameters. However, our scheme 
cannot provide a natural reason for vanishingly small #13 if it turns out so. We presented 
the results accommodating only the LMA solution, but the similar conclusion can be drawn 
also in the case of the LOW solution as can be inferred from our discussions. 


10 


Acknowledgment: EJC was supported by the Korea Research Foundation Grant, KRF- 
2002-015-CP0060, and DWJ by BK21 project by the Ministry of Education, Korea. 


[1] Y. Fukuda et. ai, Super-K collaboration, Phys. Rev. Lett. 81 (1998) 1562. 

[2] Q.R. Ahmad et. al, SNO collaboration, Phys. Rev. Lett. 89 (2002) 011301. 

[3] CHOOZ collaboration, M. Apollonio eta. ai, Phys. Lett. B420 (1998) 397. 

[4] Z. Maki, M. Nakagawa and S. Sakata, Prog. Theor. Phys. 28 (1962) 870. 

[5] L. Hall and M. Suzuki, Nucl. Phys. B231 (1984) 419. 

[6] E.J. Chun and S.K. Kang, Phys. Rev. D61 (2000) 075012. 

[7] E.J. Chun, D.W. Jung, S.K. Kang and J.D. Park, Phys. Rev. D66 (2002) 073003. 

[8] M. Hirsh, M.A. Diaz, W. Porod, J.C. Romao and J.W.F Valle, Phys. Rev. D62 (2000) 113008. 

[9] A. Abada, S. Davidson and M. Losada Phys. Rev. D65 (2002) 075010. 

[10] A.S. Joshipura and R.D. Vaidya, Nucl. Phys. B639 (2002) 290. 

[11] In light of recent SNO data, another possibility has been investigated by A. Abada, G. Bhat- 
tacharyya and M. Losada, Phys. Rev. D66 (2002) 071701. 

[12] F. Gabbiani, E. Gabrielli, A. Masiero, and L. Silvestrini, Nucl. Phys. B477 (1996) 321. 

[13] B. Mukhopadhyaya, S. Roy and F. Vissani, Phys. Lett. B443 (1998) 191; E.J. Chun and J.S. 
Lee, Phys. Rev. D60 (1999) 075006; S.Y. Choi et. ai, Phys. Rev. D60 (1999) 075002; W. 
Porod et. al, Phys. Rev. D63 (2001) 115004. 


[14] For the case of stau decay, see, M. Hirsch, W. Porod, J. Romao and J.W.F Valle, bcp- 

| ph/0207334 ; M.A. Diaz, R.A. Lineros and M.A. Rivera, hep-ph/0210182. 

[15] S. Davidson and M. Rosada, Phys. Rev. D65 (2002) 075025. 

[16] Y. Grossman and H.E. Haber, Phys. Rev. D59 (1999) 093008. 

[17] K. Choi, E.J. Chun and K. Hwang, Phys. Rev. D64 (2001) 033006. 


11 









FIG. 1: The tree-dominant points allowing the atmospheric and solar neutrino masses and mixing 
in terms of the two input variables, x\ = rn\ l /m 2 Hi and p\ = AB\/B. 



FIG. 2: Same as FIG. 1 with x '2 = m 2 L2 /m 2 Hl and p 2 = ABi/B. 
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FIG. 3: All the points allowing the atmospheric and solar neutrino masses and mixing in terms 
of the two induced variables, and £ 2 , controlling the tree-level mass matrix as in Eq. (Q). 
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FIG. 4: Same as in FIG. 3 with £2 and £3. 
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